Blow-up of the hyperbolic Burgers equation 

Carlos Escudero* 
Abstract 

The memory effects on microscopic kinetic systems have been sometimes modelled by means 
of the introduction of second order time derivatives in the macroscopic hydrodynamic equations. 
One prototypical example is the hyperbolic modification of the Burgers equation, that has been 
introduced to clarify the interplay of hyperbolicity and nonlinear hydrodynamic evolution. Previous 
studies suggested the finite time blow-up of this equation, and here we present a rigorous proof of 
this fact. 
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One of the major problems in nonequilibrium statistical mechanics is to derive an accu- 
rate hydrodynamic description of a kinetic process. A classical example is the derivation 
of the Euler and Navier-Stokes equations from the Boltzmann equation [l, 2|, for instance 



by means of the Chapman- Enskog expansion |3|, |4|- While the Navier-Stokes equations rep- 
resent a remarkable success of the theoretical study of fluid mechanics, it is well known 
that the spectral properties of their solution do not agree with experimental data for short 



wavelengths [5]. In order to solve this problem, one is tempted to derive generalized hy- 
drodynannc equations cont.nu.ng the Chapnran-Enskog expansion to higher orders Q. to 
obtain the so called Burnett equations, but they. have never achieved any notable success [7j . 
Another possibility was suggested by Rosenau [8( , who found that the telegraphers equation 

d tt u + d t u = d xx u (1) 

reproduces the spectrum of its microscopic counterpart, the persisting random walk, almost 
exactly. Using this fact, he claimed that the Chapman-Enskog expansion should be sub- 
stituted by a different expansion keeping space and time on equal footing. This procedure 
would preserve the hyperbolic nature of the resulting equations, and thus the nice spectral 
properties of the solution, at least in the linear regime 8|. An expansion of this type was 
carried out by Khonkin [9(, who found new equations for the momentum and energy fluxes 
which, in contrast to the classical Navier-Stokes and Fourier laws, depend on the first time 
derivative of these fluxes. This dependence implies in turn the appearence of a term propor- 
tional to the second order time derivative of the velocity, among others, in the corresponding 
modified Navier-Stokes equation, which becomes now hyperbolic. However, it was already 
argued by Rosenau that hyperbolicity in union with nonlinear hydrodynamical evolution 
might result in the nonexistence of the solution. A similar idea was proposed later 
where hyperbolicity was introduced to take into account memory effects in the hydrody- 
namic description of the flow, and to get rid of the infinite speed of signal propagation. In 
order to understand the interplay between hyperbolicity and nonlinear convection, a hyper- 
bolic modification of the Burgers equation (see Eq. (jSJ) below) was studied in Ref. [lO| by 
means of linear and numerical analyses. One of the conclusions of this work is that this 
equation has blowing up solutions under certain circumstances. Surprisingly, there is not 
much rigorous work done with respect to this equation. To our knowledge, there is only 
one result proving global existence in time of the solution provided the initial data is small 



id, 
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enough [ill ]. It is thus our goal to extend this result and prove blow-up of the solution for 
large initial data. But before starting with the analysis, let us discuss the validity of the 
hyperbolic Burgers equation as a physical model. It is important to note that the hyperbolic 
hydrodynamic equations derived from the Boltzmann equation [9| and obtained including 
memory effects [10(] are far more complex than the equation under study. Furthermore, 
there is, to our knowledge, no direct application of the hyperbolic Burgers equation to a 
physical problem. It is thus to be interpreted as a model equation resembling some of the 
characteristic properties of the full models, which will hopefully help us to understand the 
dynamical properties of the more difficult equations. 

In this article we are concerned with the Cauchy problem for the hyperbolic Burgers 
equation 

fid a v + d t v + vd x v = ud xx v, (2) 
v(x, 0) = v (x), (3) 
v t (x,0) = v 1 (x). (4) 

where v represents the one-dimensional velocity and fi > and v > stand for the inertia 
and the viscosity respectively. From now on we will assume that the initial conditions are 
compactly supported in the interval J Q = [— L, L], for some < L < oo. We will use the 
notation c = a/ vj fi, because, as we will see, this quantity will play the role of the speed of 
sound in Eq. (j2J), that is, the maximum velocity of propagation of disturbances. The proof 
presents some similarities to the corresponding one of nonexistence of the solution to the 
compressible Euler equations jl^ . but, interestingly, the singular behaviour developed by 
the hyperbolic Burgers flow is not of the shock-wave type. Actually, what diverges in this 
case is the velocity itself, as we will show. 

To prove the finite time blow-up of Eq.(j2D we need first some result concerning the finite 
velocity of propagation of disturbances. 

Proposition. Let C(x,t) be the cone 

C(x,t) = {(x,t):\x c -x\<c(t c -t)}, 0<t<t c , (5) 

for some x c G WL, t c E M + . If v is a classical solution of Eq.fijQ) in C(x,t) and v , v x = in 
C(x, 0), then v = m C(x, t),0<t<t c . 
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Proof. Let us rewrite Eg. (12]) in the form 



d tt v + -d t v + -vd x v = -d xx v, (6) 



and define the energy functional 



E{t) = \ l^ + c 2 v 2 x )dx. (7) 



2 

The first derivative of the energy functional yields 

E\t) = / (v t v tt + c 2 v x v xt )dx - %o 2 + c 2 v 2 x } ac - - ^[v 2 + c 2 v 2 ] dc +, (8) 

where dC +( -~^ denotes the right (left) boundary of C(x, t) for a fixed time, and an expression 
between square brackets with this subindex indicates that it is evaluated at this point. We 
can further proceed integrating by parts 

f c c 

E\t) = J v t (v tt -c 2 v xx )dx+c 2 [v x v t } dc +-c 2 [v x v t } dc ~--[v 2 t +c^ (9) 

and using Cauchy inequality in the boundary terms 



c 2 [v x v t } dc + - c 2 [v x v t ]ac- < | {[v 2 + c 2 v 2 x ] dc + + [v 2 + c 2 v 2 x ] dc -) , (10) 



in order to get 



E'(t) < [ v t (v tt - c 2 v xx )dx = / v 2 dx / v t vv x dx 

Jc V Jc V Jc 

< 1 / \v t VV x \dx < ML-^L-iC)) f {v , + c2vl)d ^ (n) 

fl J Q ZfiC 



J c 

where we have used Cauchy inequality together with the regularity of classical solutions to 
Eq. (j2|). We can integrate this last expression to obtain the inequality 

E(t) < exp ( M^^m t ) E(0), (12) 



which, taking into account the fact that E(0) = 0, yields the identities Vt, v x = in C(x,t), 
and thus, in addition to the initial conditions, the desired result v = inside the cone 
C(x,t). □ 

Remark. Note that this result also holds for v being a bounded weak solution of Eq. (0|) 
in C(x, t). 
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Let us now define Jclas 

J{t) = {x e R : \x\ < L + ct}. (13) 

As a consequence of the proposition we know that if our initial conditions are compactly 

supported in the interval Iq, then a classical solution to our problem will be compactly 

supported in J(t). The proof of our main result will rely on constructing a differential 
inequality for the quantity 

F(t) = / xvdx, (14) 



which is a measure of the expansion of the velocity. We know that this quantity is convergent 
as long as the solution is regular enough, by direct application of the proposition. 

Theorem 1. Suppose v is a classical solution of Eq.$$ for < t < T. If the initial 
conditions fulfill 

16 

F(0) > — cL(L + 6cfi), (15a) 
3 

64 

F'(0) > — c 2 (L + 6c/i), (15b) 

then the life span of the classical solution is finite. 

Proof. Consider the effect of the time evolution on F(t): 

I^F"{t) + F'(t) = I x(fid u v + d t v)dx = f x{—vd x v + vd xx v)dx — - j v 2 dx, (16) 

where we have used integration by parts and the compact support property of the classical 
solutions to Eq.(j2j). Now we derive 

xvdx 1 = I / xvdx I < / x 2 dx / v 2 dx = 
) \Ji(t) ) Ji(t) Ji(t) 

-{L + ctf [ v 2 dx = -(L + ct) 3 [v 2 dx, (17) 
3 Ji(t) 3 J K 

where we have used again the fact that the solution is compactly supported together with 
Schwartz inequality. By using these two relations we obtain the following differential in- 
equality 

^F"{t) + F'{t)>h y L + ct)- 3 F{t) 2 . (18) 



We can now adapt the argument in Ref. [13j to prove the blow-up of F(t). Consider the 
auxiliary initial value problem 

G'(t) = e(ct + L)- 3 G(t) 3/2 , (19) 

G(0) = F(0), (20) 



where e > is to be set later. It can be solved to yield 

G(t)- 1/2 = G(0)-^ + ^ [(* + L/c)-" - {L/c)-'] . (21) 

This solution blows up in a finite time T* provided the initial condition fulfills 

G(0) > (22) 

and in this case it also increases strictly and monotically from the initial condition to infinity. 
Apply now a time derivative on Eq. (fT9l) to find 

G"(t) = ^-e 2 (ct + L)- 6 G(t) 2 - 3ec(ct + L)^G(t) 3/2 < |e 2 (ct + L)- 6 G(t) 2 . (23) 

Summing this last equation multiplied by \i and Eq. (fT9]) we get 

fiG"(t) + G'(t) < ^e 2 (ct + L)- Q G{tf + e(ct + L)- 3 G(t) 3/2 , (24) 



that yields in turn the inequality 



nG"{t) + G'(t) < 



eG(0) " /2 + ^ e2 



(ct + L)- 3 G(t) 2 . (25) 



Now choose e small enough so that the relations 

3_/i 2 <• 3 
2L 3 - 4 

G'(0) = eL- 3 G(0) 3/2 < F'(0), (27) 



eG (°)" 1/2 + if^ 2 - z> (26) 



hold. So we have finally reduced our problem to the inequality 

fiG"(t) + G'{t) <\{ct + L)- 3 G(t) 2 , (28) 

with initial conditions fulfilling G(0) = F(0) and G'(0) < F'(0). Now, it only rests to prove 
that F(t) > G(t). The fact that F'(0) > G'(0) implies that there exists some time t > 
such that F'(t) > G'(t) if < t < t Q . Let us denote by ti the supremum of all the values of 
t . Suppose tx < T\ then F{t) > G(t) for < t < t h and F'fa) = G'(ti). Thus we have 
F'(t) - G'{t) > implying that F(t) - G(t) > F{0) - G(0) = for < t < t v This, in turn, 
implies that F(ti) > G(ti), while the difference between Eqs. ( ITBl and ( 1281 yields 

li [F"(t) - G"{t)\ + F'it) - G\t) > l(ct + L)- 3 [F(t) 2 - G(t) 2 } > 0, (29) 
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for < t < t\. This inequality is equivalent to 

j t (e*/" [F'(t) - G'(t)]) > 0, (30) 

which can be straightforwardly integrated in order to get 

e *i//i [F'(ti) - G'{t x )] > F'(0) - G'(0), (31) 

or, what is the same, F'(ti) > G'(ti). So now we can argue by contradiction that ti > T* 
to conclude the proof. □ 

Remark. Note that the requirements on the initial conditions Eqs. U5a\) and U5b\) come 
from Eqs. fifty) , ( ESP and fiFfy in the proof. 

These arguments suggest that the velocity is blowing up in finite time, but they do not 
constitute a proof of blow-up, just of nonexistence of the classical solution. This is because 
the maximal time of existence of the classical solution might be strictly smaller than the 
blow-up time, as happens in certain situations fl3 |. To be sure that the velocity increases 
unboundedly in finite time we need an additional theorem stating that the solution can be 
continued as long as we have controlled the L°° norm of v. For this we can assume the 
existence of a weak solution to Eq.p]); interpreting all the derivatives in the weak sense 
leaves us with the requirement v G L 1 fl L 2 ([0,T] x K) for the weak solution to exist. This 
fact, in addition to the proposition, implies that the nonexistence of the weak solution 
is necessarily due to a divergence of the velocity. Before starting, let us define the set 
n = (0,T) x U{t) C M 2 , where U{t) is open, bounded, J{t) C U{t) C R, and such that 
Q has a C 1 boundary. Henceforth we will assume that the initial conditions are compactly 
supported in the interval 1$ and they fulfill v$ G H 3 (M.) and v\ G H 2 (M.). 

Theorem 2. Suppose there exists a weak solution v to the Cauchy problem for Eq.ftty) 
for < t < T < oo. // this solution fulfills \ \v ||z,°°(n) < °°; then it is a classical solution. 

Proof. First of all, we know that a bounded weak solution of Eq.([2]) has compact support 
as a consequence of the proposition. Then, it is very easy to see that 

IMb(n)= / / v 2 dxdt<2\\v\\*)(^T 2 + LT). (32) 

JO JR XZ J 

We can now evaluate the quantity 

E x {t) = \ [ (v? + c 2 vl)dx, (33) 



following a parallel reasoning to that of the proposition, to get 



This allows us to conclude 

f -T 



[v 2 + c 2 v 2 x )dxdt < E(0)- 



'0 JR 

Consider now the functional 



MU°°(fi) 



exp — T 

/iC 



(34) 



(35) 



(36) 



Its first derivative is 



^(0 = / {vttVttt + c 4 v xx v xxt )dx = / (v tt v ttt - c 4 v xxx v xt )dx, (37) 

JR JR 

after derivation by parts. Now we will proceed to evaluate the second integral in the right 
hand side of the equation above 



^ / ^xxx^xtdx 



Vxt 



1 1 / \ 

Vttx + ~V xt + -{VV X ) X 



We have 



h 



1 d 
Ydt 



v x tv xtt dx = ™ / vl t , 



Kt > 0, 



and 



1 f 1 2 \ 1 d / q , ||'^'||L oo (n) 1 1 ■■> ■> , 
h = - I [v xt v x + v xt vv xx )dx < — — I v x dx H — — / (v xt + c v xx )dx. 



The first integral in the right hand side of Eq. (137j) can be estimated as follows 



VuVatdx = / v tt 



—Vu K)( + c 2 v. 



xxt 



dx = J\ + J 2 + J3. 



We have 



// 



%cfe < 0, 



and 



•/;! : I '',/''.,-,/'/■'• • ~r 2 / v xt v xtt dx = --^^ I ''^(h: 



(38) 

(39) 
(40) 

(41) 

(42) 
(43) 



h = -- / {vttVtV x + v u vv x t)dx < — - (n) / (^ + c 2 v 2 xt )dx - -~ / (v 2 ) t v x dx, (44) 



(45) 



Now we can integrate Eq. (!37|) with respect to time and use the above expressions to obtain 

[y 2 xt + c 2 v 2 xx )dxdt 



( v tt + cv xx + 2c v xt )dx 
3/i 



c|M|z,«(n) rU 



< 





Ti 


v\dx 


+ 







1 






r 


+ 27f 


Jo 



2/i 

\ v \ \l°°(q) 



rTi 



V c Jo 



v 2 v xt dxdt — 



(4 + c 2 v 2 xt )dxdt 

i 

r 1 Ti 



2/i _ 



v t v x dx 



(46) 



after integration by parts of the last term. We still have to estimate the following terms 



/ / v 2 v xt dxdt = - / 
'o Jr & Jo 



(v t ) x dxdt = 0, 



(47) 



due to the compact support property, 

77— / vldx = -77- / vv x v xx dx < - \ L (0) 
Jr 3/i J R 3/i 



for e > 0, and 
1 

2jl 



v\v x dx = / vv t v xt dx < - ^ L — 

^ Jr V 



+ ) d.r. 



1 2 ^ 2 1 7 
Ys v t + ~ 2 v xt ) dx. 



(48) 



(49) 



for 5 > 0. Now we can choose e and 5 small enough, say 



e < 



3/ic 2 



5 < 



2/ic 2 



oil 11 11 I • (50) 

^||f||L°°(n) 1 1^| \L ao (Q) 

and noting that for an arbitrary Ti, such that < Ti < T, we have | | |i°°(o,T!) < dor 
some suitable constant C, we can rearrange inequality (j4"6l) in the following way 

{vl + c 2 t; 2 t + c 4 v 2 xx )dx] (Ti) <C 1 + C 2 [ 1 [ {v 2 t + c 2 ^ t + c 4 v 2 xx )dxdt, (51) 
for some suitable positive constants Ci and C2, and recall Gronwall inequality to find 



/ 2 1 22 1 4 2 \ j 

Finally, integrating with respect to time yields 



{n) <Ciexp(C 2 Ti). 



(52) 



/ T / (4 + c 2 t; 2 t + cS 2 J^t < (<7i/<7 2 )[exp(C7 2 T) - 1]. (53) 



Now we know that the Sobolev norm 

r-T 



Ml# 2 (n) 



L<>0 



{(% + cv xt + c v xx )n + (v t +c v x )fi + v }dxdt 



1/2 



(54) 
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is finite, so we can use a Sobolev embedding to obtain 

max(| IubIIl^q), |H|£°o ( fi)) < oo. 

This inequality will prove itself very useful in order to estimate the functional 

Es{t) = \ / K t + c & v 2 xxx )dx. 
z Jr 

Its first derivative reads 

E 3(t) = / (vutVtttt + c 6 v xxx v xxxt )dx = / (vtttVtttt - c 6 v xxxx v xxt )dx. 



Let us start with the first integral in the right hand side of this equation 



vtttvttttdx = — I v 2 tt dx - - / v t tt{vv x )udx + c 2 I v m v xxtt = K x + K 2 + K 3 . 
M7r M 



We see that K x < 



K 2 = -- [ [v ttt (v tt v x + 2v t v xt + vv xtt )\ < Vx L 2 {Q) [ (fi 2 v 2 tt + v 2 t )dx + 
^ Jr l[i 



V c J R 2/ic 

and 

1 d 



K 3 = - v xttt v xtt dx = --— v xtt dx. 
Jr * at Jr 

The second integral in the right hand side of Eq.( [571) may be evaluated as follows 



c I v xxx v xxx tdx c I v xx iV xxxx dx 
Jr Jr. 

/ {v xxt [v xx tt + ^~ X v xxt + n~ l {vv x ) xx \}dx = Li + L 2 + L 3 
Jr 



We see that L 2 < 0, 

1 d 
Ydt 

and 



L\ — — --77 / v 2 t dx, 



t 1 f /o \ 1 ^ 3 Pi L»(0) [ / 2 2 2 \ 7 

L 3 = / v xxt (3v x v xx + vv xxx )dx < — / (/i + v xx )dx 

H Jr Jr 

+ V L °° in) [ (v 2 xxt + c 2 v 2 xxx )dx. 



10 



Collecting these results, we can transform Eq. ljSTj) into the inequality 
d 



dt 



%t + c v xtt + c v xxt + c v xxx )dx < D 1 / (v tt + cv xt + c v xx )dx + 



I 2 , 2 2 i 42 i 6 2 \ j 



(64) 



for suitable positive constants D\, D 2 . Now, integrating this equation with respect to time 
we obtain 



I i i 22 i 42 , H 2 \ t 

(%t + c V xu + c v xxt + c v xxx )dx 

T 2 



(T 2 )< 



Di + D 2 



v t tt + c v xtt + c v xxt + c v xxx )dxdt, 



(65) 



for < T2 < T, where we have used the boundedness of the Sobolev norm 1)541) . Let us 
recall once again Gronwall inequality to get 



/ 2 1 2 2 ,42 1 e 2 \ j 

i V ttt + C V xtt + C V xxt + C V xxx)dx 



(T 2 )<D ie x V (D 2 T 2 ). 



(66) 



Integrating this expression with respect to T 2 in the interval [0, T] yields 

/ / « + c 2 v 2 xtt + c 4 v 2 xxt + c% 2 xxx )dxdt < 00, (67) 

JO JR 

and now, if we put this result together with the boundedness of ||u||H 2 (fi) we obtain the 
boundedness of the norm ||u||ir 3 (n) which implies, after invoking the corresponding Sobolev 
embedding, the desired result v G C 2 {VL). □ 

Corollary. Consider the Cauchy problem for Eq.^Q) with initial conditions fulfilling 
the same properties as in Theorems 1 and 2. Then there is a finite time t m such that 
rim f _* tm |M|l°°(b 



00. 



In summary, we have proven two results that were previously conjectured in the physics 
literature concerning the hyperbolic modification of the Burgers equation. The first one es- 
tablishes the finite speed of propagation of disturbances in this type of flow, and provides us 
with an estimation of the velocity of sound within it, c = \Jvj y«. The second one shows that 
the solution to the hyperbolic Burgers equation blows up in finite time, provided that the 
initial condition is large enough, a result that complements the partial regularity theorems 
proven by Rudjak and Smagulov In contrast to what happens in the case of an inviscid 
compressible fluid [12J , the flow does not develop a shock- wave, but the velocity itself diverges 
in finite time. Physically, we have shown what effects have in the macroscopic description 
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including memory in the corresponding microscopic model. Although the undesirable prop- 
erty of infinite speed of signal propagation is lost, now we find that the nonexistence of the 
solution complicates the use of hydrodynamic equations to describe the coarse-grained dy- 
namics of the microscopic process. Furthermore, the fluid might self-accelarate till reaching 
an infinite velocity, a fact that is clearly nonphysical, and that rules out the possibility of 
using the hyperbolic Burgers equation for large initial conditions. We hope that the analysis 
of this model will facilitate the assessment of the more complex hydrodynamic equations 
that appear as asymptotic approximations of the Boltzmann equation when a regularized 
version of the Chapman-Enskog expansion is performed. 
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problem. This work has been partially supported by the Ministerio de Education y Ciencia 
(Spain) through Projects No. EX2005-0976 and FIS2005-01729. 
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